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Abstract. We give a computationally efficient method for constructing the linear dif- 
ferential operator with polynomial coefficients whose space of holomorphic solutions is 
spanned by all the branches of a function defined by a generic algebraic curve. The 
proposed method does not require solving the algebraic equation and can be applied in 
the case when its Galois group is not solvable. 



1. Introduction 

To find relations satisfied by a given special function is a difficult and im- 
portant problem in the theory of special functions of mathematical physics. 
The relations in question can involve derivatives, integrals, finite differences 
etc. Knowing a global relation for a special function that is defined locally 
(e.g. by means of a series converging in a neighbourhood of a point) allows 
one to deduce global properties of that function. From this point of view, 
linear differential equations with polynomial coefficients are of particular 
interest. One of the reasons for this is the difficult problem of computing 
the analytic continuation along a given path of a locally defined special 
function. By identifying such a function with a solution to a system of lin- 
ear differential equations with polynomial coefficients which does not have 
any "extra solutions" one can use standard techniques for investigating the 
analytic continuation of the function under study (see, for instance, |H]). 
Here by "extra solutions" we mean the solutions which are not branches 
of the function under study, that is, which cannot be obtained from it by 
means of analytic continuation. Observe that every germ of a (multivalued) 
analytic function satisfies a relation with entire (in particular, polynomial) 
coefficients provided that this relation is valid for one of its germs in a 
neighbourhood of some fixed nonsingular point. 

The culmination of this approach is the Wilf-Zeilberger algorithmic proof 
theory (see [IE] and [12]) based on holonomic systems of equations. In the 
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present paper we thoroughly investigate the special case when the func- 
tion under study is algebraic and the holonomic system consists of a single 
ordinary linear homogeneous differential equation with polynomial coeffi- 
cients. Despite all simplicity, this setup leads to formidable computational 
challenges. 

The 21st problem in the Hilbert list was solved in 1989 by A. A. Bolibrukh 
who proved that it is in general not possible to construct a linear fuchsian 
system of differential equations with a prescribed monodromy group (see 
PQ). However, the problem of effective computation of a system of differ- 
ential equations (and, in particular, of a single differential equation) with a 
prescribed branching of solutions (whenever this is possible) remains open 
and is in the focus of intensive research, see [3], [6], and The com- 
puter algebra system Magma has a built-in command DifferentialOperator 
for finding such operators (see [5]). In the paper [11] a method for com- 
puting annihilating operators for a class of algebraic functions was devel- 
oped. However, due to computational difficulties and software limitations, 
it could only be used for computing annihilating operators for algebraic 
functions defined by equations with solvable Galois groups. In the present 
paper we describe an algorithm which allows one to compute annihilating 
operators for an essentially larger class of algebraic functions and does not 
require the solvability of the Galois group (see Examples [H (4), E| dU] and 
O). 

It is well-known that an ordinary linear differential equation with a pre- 
scribed solution space can be found by means of the wronskian of a basis of 
this space. However, from the computational point of view, the wronskian- 
based representation of the differential equation for an analytic function 
(which is, in general, defined only locally) is merely a nonconstructive ex- 
istence theorem. There are three main reasons for this. First, to form 
the wronskian, one needs to choose a basis in the space of germs of the 
given function at a nonsingular point. This requires computing the analytic 
continuation of the given function along any path which is, in general, a dif- 
ficult problem. Secondly, to evaluate a determinant containing high-order 
derivatives of a given special function is a task of a great computational 
complexity. Finally, extracting the polynomial coefficients of the desired 
differential operator out of the obtained combination of algebraic functions 
requires the full use of modern methods of computer algebra. For instance, 
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to compute the differential operator for the roots of the generic monic cu- 
bic by means of the wronskian is already a challenge (see example in §5 in 
[IT]). In the general case, the wronskian-based construction is not suitable 
for computation since no effective means of simplifying expressions which 
contain high-order derivatives of special functions are presently known. 

The present paper provides an algorithm for constructing the optimal 
(that is, of the smallest possible order) linear homogeneous differential 
equation with polynomial coefficients for a univariate algebraic function 
y = y(x) implicitly defined by the equation 

(1) y m + a iy m i + ... + a n y m «+x = 0. 

This algorithm allows one to overcome the difficulties listed above. That is, 
it does not require to solve the problems of analytic continuation, evalua- 
tion of determinants and simplification of expressions involving derivatives 
of algebraic functions. The proposed method is a development of the ideas 
of the work [11] . It allows one to reduce the problem of computing the 
annihilating operator for an algebraic function to the problem of finding a 
basis in the syzygy module of an ideal in the ring of multivariate polynomi- 
als. The presented algorithm differs from other methods (both recent and 
classical, see pf], 0, [S], and P]) in its primary field of application (it deals 
with generic algebraic equations), in the underlying concept (holonomic 
systems of partial differential equations and noncommutative elimination) 
and the complexity of differential operators that it can efficiently produce. 
The capabilities and limitations of the proposed algorithm are summarized 
in Table 1. 

The authors are thankful to Professor D. Zeilberger for helpful expla- 
nations giving insight into holonomic systems approach and to Professor 
M. Singer for comments on Galois theory. 



2. Annihilating operators for solutions of holonomic 
systems of differential equations 

In what follows we will denote by V n the Weyl algebra of differential 
operators with polynomial coefficients in n variables x = (x±, . . . , x n ) G C". 
This algebra is generated by the operators x\, . . . , x n: . . . , satisfying 
the relations ^ o Xj — Xj o = 5^. Here "o" denotes the composition of 
differential operators. The Weyl algebra is simple (see Chapter 1 in [2]). 
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When speaking about ideals in the Weyl algebra, we will always mean its 
left ideals. 

The following basic statement is well-known but not easy to find in the 
literature in the following explicit form. It can be deduced from Theo- 
rem 2 in [TTJ. It also follows from Theorem 1.4.12, Proposition 1.4.9, and 
Lemma 2.2.3 in [15]. 

Proposition 1. For any holonomic left ideal I C T> n and any i G {1, . . . , n} 

there exists a nonzero operator Pi G /, all of whose derivatives are with re- 
spect to the variable X{, that is, an operator of the form 



, dx] 



The following statement is a consequence of the results in |TD] and 
It shows that algebraic functions defined by generic algebraic curves are 
annihilated by holonomic ideals in T> n . 

Theorem 2. Any germ of the algebraic function y(xo, X\, . . . ,x n ) implicitly 
defined by the relation 

(2) x n y n + x n . x y n ~ l + . . . + x x y + x = 0, 

satisfies the holonomic system of differential equations 



(3) 



J* S = J* I . whenever i + j = k + /, 

n n 

E*»i^ = -!/ and E^aj = 0. 



i=0 i=0 

Conversely, any holomorphic solution of (0j defined locally in a neighbour- 
hood of a nonsingular point is a linear combination of germs of the function 
y(xo, xx, . . . , x n ) at this point. 

The system of differential equations (ED is a special instance of the Gelfand- 
Kapranov-Zelevinsky hypergeometric system introduced in ||10j |. Its "de- 
homogenized" version for an algebraic curve with affine parameters was 
investigated by Mellin in [T3] . 

Recall that the Nilsson class comprises (multi-valued) analytic functions 
of several complex variables which have finite determination and moderate 
growth in arbitrary neighborhood of any of their singularities (see §4.1.12 
in [3]). Here by the determination of a multi-valued analytic function we 
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mean the number of its linearly independent germs in a neighborhood of a 
generic point in its domain of definition. The determination of an analytic 
function of one complex variable which lies in the Nilsson class and has 
finitely many singularities in C coincides with the smallest possible order 
of an ordinary linear homogeneous differential equation with polynomial 
coefficients satisfied by this function. 

Theorem E] together with Proposition Q] imply the existence of a lin- 
ear differential operator with polynomial coefficients whose space of local 
holomorphic solutions at a nonsingular point is spanned by the roots of the 
generic algebraic equation (ED and all of whose derivatives are with respect 
to Xq. This operator is defined uniquely up to a sign. We will say that this 
operator is optimal for the given generic algebraic curve. 

Example 3. Consider the algebraic function y(xo,xi,X2) defined as the 
solution to the quadratic equation X2I/ 2 + x\y + xq = 0. By Theorem [21 any 
of its branches lies in the kernel of any operator in the ideal J with the 
generators 

d 2 d 2 d d d d d 
A = - — = xi- \-2x 2 ^ hi, 6 = x - hxi- hx 2 — • 

OXqOX2 UX\ OX\ OX2 OXq OX\ OX2 

Since the determination of the function y{xQ,x\,X2) equals 2, Proposition 
CD yields the existence of a second order differential operator P G J, all of 
whose derivatives are with respect to xq. Using the notation 9i = X{-^, we 
can write the expansion of this operator with respect to the basis of J in 
the form 

P = Xqx\x2A — {{x\ — 2^0X2)^0 + XqX29\)B + {{x\ — 4X()X2)Oo + 2xqX20\)C 

_ iX0X2) f 2 _ 2X2 JL) . 

OXq OXq J 

Of course, this optimal differential operator is only a monomial multiple of 
the wronskian of the roots of the initial algebraic equation. 

In the next section we describe the algorithm for computing the opti- 
mal annihilating operator for an arbitrary algebraic function satisfying an 
equation of the form (JTj). This will, in particular, perform the noncom- 
mutative elimination of all the derivatives except for -M^- in the holonomic 
ideal ([3]) by means of methods of commutative algebra only. 




6 



V.A. KRASIKOV AND T.M. SADYKOV 



3. Computing the annihilating operator for a given 



ALGEBRAIC FUNCTION 



We begin by computing the determinations of some elementary functions 
and the corresponding differential equations. 

Example 4. (1) Any rational function / = p(x)/q(x), where p(x), q(x) G 
C[x], has determination 1 and satisfies the first-order differential equation 

pqf = {p'q-pq')f- 

(2) The function / = x a also has determination 1 for any a G C since 
its analytic continuation e 2nta f around the only finite singularity x = is 
proportional to /. It satisfies the differential equation xf = af. 

(3) The function / = lux has determination 2, since its analytic continu- 
ation along any path can be written in the form In a: + 2irki, k G Z. Thus 
any germ of / at a nonsingular point lies in the two-dimensional linear 
space with the basis {l,lnx}. The second-order differential equation with 
polynomial coefficients satisfied by / has the form xf" + f = 0. 

(4) The algebraic function y = y(x) implicitly defined by the relation 
y h + ay + x = has determination 4 (see Theorem below) and satisfies 
the differential equation (256a 5 + 3125x A )y^ + 31250x 3 y^ + 73125xV 2 ) + 
31875x2/' - 1155?/ = 0. 

(5) Finally, the function 1/lnx has infinite determination since its germs 
{l/(lnx + 27r/cz)} fcgZ are linearly independent. This implies, in particu- 
lar, that this function does not satisfy any linear homogeneous differential 
equation with polynomial coefficients. 

In the present section we describe an algorithm for computing the optimal 
annihilating operator for the roots of a generic algebraic equation with 
symbolic coefficients, that is, an equation of the form (00). The roots of 
the equation aoy m + a\y mi + a2y m2 + . . . + a n y mn + a n+ \ = (regarded as 
functions of a = (ao, . . . , a n +i)) satisfy the holonomic 74-hypergeometric 
system with the vector of parameters (0,-1) (see [IS]), where 



Namely, it is the left ideal in the Weyl algebra C[ao, . . . , a n +i 3 do, ... , d n +\\ 
generated by 

the toric operators d u — d v for u, v G N n+2 with A • u = A • v, 




1 1 

rru 
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n+1 n 

and the Euler operators djdj and mao + rrijajdj + 1. 

3=0 j=l 

Thus by Proposition Q], there always exists a linear differential operator 
with polynomial coefficients in clq, . . . , a n +i, and all of whose derivatives 
are with respect to a n+ \. Setting ao = 1 and a n+ \ = x we obtain the 
annihilating operator for the solutions of (ED). Using noncommutative elim- 
ination theory, one can compute this operator in a way similar to that in 
Example El In the special case of a trinomial equation (that is, for n = 1) 
the desired operator is a right factor of the Mellin differential operator 
found in P~5| . 

The following theorem gives the order of the annihilating operator. 

Theorem 5. (E. Cattani, C.D'Andrea, A. Dickenstein [7]) The number 
of linearly independent (over the field of complex numbers) germs of the 
solutions to the equation (QP at a generic point x G C and for generic 
values of the parameters (a±, . . . , a n ) G C n is given by 

if GCD(ra,mi, . . .,m n ) = 1, 
K{m, mi, ... , m n ) - < m f GCD ( ) L 

k GCD(m,mi,...,m„)' ^ v j xj in/ 

Here [ } denotes the integer part of a real number. 

The following theorem is the foundation of our algorithm for computing 
optimal annihilating operators. 

Theorem 6. Let s, = ), i = 1, . . . , m be the roots of the 

algebraic equation (QP and denote P(t) = n^i^ — s «)- -^ or e,yer ?/ k = 
1, . . . , m we define the ideal Ik in the polynomial ring with m + n + 1 
variables C[si, . . . , s TO , ai, . . . , a n , x] to be 



■D'(^-l)!| \ \(s k - Si) 2 ™' 1 | res-^, £ = 1, . . . , i?(m, m 1? . . . , m r 

t=S K P(t) 




The vector of polynomial coefficients of the optimal annihilating operator 
for the algebraic function defined by (Ep lies in the following syzygy module 
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of the quotient of the ideal Ik with respect to the Vieta relations: 

Syz(7 fc / ( S m - m .(si, . . . , s m ) - (-l) m ~ m i a h j = 1, . . . ,n; 
(4) S m -k{si, ■ ■ . ,s m ), for k £ {0, mi, m 2 , . . . , m n }; 

S m (si,...,s m )-(-l) m ar)). 

Here <Sj(si, • • • , s m ) 25 £/ie elementary symmetric polynomial of order j in 
the variables (si, . . . , s m ). In the sequel we will denote the ideal generated 
by the Vieta relations by V. 

Proof. Let x G C be a point outside of the zero locus of the discriminant 
of the left-hand side in (CD). Let (a±, . . . ,a n ) G C n be a generic vector of 
parameters and let ai, . . . , a n ) denote the k-th branch of the solution 
y(x, ai, . . . , a n ) to the algebraic equation (Q]). Let us now denote by D the 
differential operator 

D _ 5 ™ 

dsi...ds m ' 

Using the well-known contour integral representation for a solution to a 
univariate algebraic equation (see Section 5 in fTl"] ) we conclude that the 
generators of the ideal Ik are polynomial multiples of the derivatives of the 
solution to (0Q): 

d e y k {x,ai,...,a n ) 

= (-w- l )'a* = i^M" 1 (sum) 

- ((^-l)!)— ^(* fc -ai)...[fc]...(a*-*„)y 



((s fc -ai)...[A:]...( Sfc - am ))' 

= _ v- (e+i 1 -iy....[k]...{e+i n -i)\ 

n ^ ((£-l)\) m - 2 hl..[k]...i m \ {s k -siY+n...[k}...(s k -s n y+^ ■ 

ii+...[k\...+i n =£-l 

This shows that the generators of the ideal Ik are indeed elements of the 
ring C[si, . . . , s m , a\ 1 . . . , a n ,x]. By Theorem \5\ the determination of the 
solution to (ED) equals R(m, mi, . . . , m n ). Thus by Theorem [2] and Proposi- 
tion Q] there exists a linear differential operator with polynomial coefficients 
(in x, a±, . . . , a n ) all of whose derivatives are with respect to x and whose 
space of holomorphic solutions at a generic point is spanned by the branches 
of y(x, ai, . . . , a n ). For the sake of computational efficiency we factor out 
the Vieta relations. This increases the number of variables involved in the 
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generators of the ideal but decreases their degrees. The desired differential 
operator is a relation between the derivatives of y{x, a±, . . . , a n ) with poly- 
nomial (and thus single-valued) coefficients. By the conservation principle 
for analytic functions the same relation must be satisfied by any of the 
germs of y(x,ai, . . . , a n ) at a nonsingular point. Thus the coefficients of 
this relation lie in the syzygy module (Hj). □ 



Observe that the elements of the syzygy module (J3| are polynomial vec- 
tors whose entries in general depend on all of the variables si,...,s m , 
<2i, . . . , a n , x. The proof of Theorem El implies that there exists an element 
of (HI) whose entries only depend on a\, . . . , a n , x. It can be found by means 
of the following algorithm. 

Algorithm 7. The actual computation of annihilating operators for alge- 
braic functions was organized as follows: 

1. Compute the basis of the ideal I\ defined in Theorem El 

2. Using the lexicographic order of the variables s\, . . . , s m compute the 
Grobner basis of the ideal V defined by the Vieta relations (as defined in 
Theorem E])- 

3. Perform polynomial reduction of the generators of the ideal I\ by 
means of the Grobner basis of the ideal V. That is, at this step, we use 
the Vieta relations as much as possible in order to simplify the generators 
ofV. 

4. Factorize the obtained family of polynomials. The result has a very 
specific structure: it is a family of polynomials in C[si, . . . , s m , a\, . . . , 
a n , x] whose elements are symmetric with respect to S2, . . . , s m . Using the 
Grobner basis of the ideal V, reduce them to polynomials in C[si, ai, . . . , 
a n ,x]. Let us denote this family of polynomials by Hi, . . . , 7Z m . 

5. Any C[ai, . . . , a n , x] -linear relation for the family of polynomials TZi, 
. . . , lZ m transforms into a linear system of algebraic equations over the 
field of rational functions in the variables ai, . . . , a n , x. Proposition Q] and 
Theorem yield the existence of an at least one-dimensional C- vector space 
of solutions to this system of linear equations. Finding a basis in this 
space and clearing the denominators, we obtain the desired polynomial 
coefficients of the optimal annihilating operator for the initial algebraic 
function. □ 
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Example 8. The linear space spanned by the roots of the algebraic equa- 
tion 

(5) y 5 + 2y A -3y 3 + y 2 + 5y + x = 

(in a neighbourhood of a point where the discriminant of this equation 
does not vanish) coincides with the linear space of holomorphic solutions 
to the differential equation 

(-43728190560 + 795819153 x - 53446888 a; 2 + 56028 x 3 ) 
(-1585575 + 71982 x + 281583 x 2 + 81342 x 3 + 3125 x 4 ) y^ + 
15(-650327879439783 - 5747872136026563 x - 2400588229818366 x 2 - 
91559102743545 x 3 - 304019551433 x 4 - 131338505212 x 5 + 128397500 x 6 ) 
60(-1821690090417321 - 1560609625036728 x - 98711280942848 x 2 + 
721492325057 x 3 - 103787727624 x 4 +91045500x 5 ) + 
180(-282046871305467 - 38794189010031 x + 478890241959 x 2 - 
28458003540 x 3 + 21944300 x 4 ) + 

720(- 1756652589603 + 23053844253 x - 812236372 x 2 + 522928 x 3 ) y' = 0. 

Observe that the Galois group of the algebraic equation (Ej) is not solvable 
and hence its roots can only be expressed in terms of special functions (e.g., 
functions of hypergeometric type, see [16]). Despite this fact, Algorithm 
[7] allows one to compute the annihilating operator for the roots of the 
equation ((3) using only the methods of commutative algebra. 

The following example provides a fundamental system of solutions to a 
fifth-order linear differential equation with polynomial coefficients. 

Example 9. For any a G C* a basis in the space of holomorphic solutions 
to the differential equation 

((256/5) a 5 x 3 + 625 x 4 )y {5) + (384 aV + 6875 x 3 )y {A) 

+ (624 a 5 x + 19500 x 2 )y {3) + (168 a 5 + 14100 x)y {2) + 1344?/ = 

in a neighbourhood of a generic point x G C is given by the roots of the 
algebraic equation y 5 + ay 4 + x = 0. 

4. Software, hardware and examples 

Most of the examples in this paper were computed by means of a Mathe- 
matica 6.0 package developed by the authors and run on an Intel Core(TM) 
Duo CPU clocked at 2.00GHz. The bottleneck of the algorithm is com- 
puting the syzygy module of an ideal in a ring of polynomials in several 
variables. In some cases, we have used Singular for this. 
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Example 10. Generic quintic. One of the goals of the research presented 
in this paper was to provide a computationally efficient extension of the 
results of Section 5 in [TT] beyond the class of algebraic equations with 
elementary solutions. In this example, we demonstrate the efficiency of 
the described approach by means of the generic monic quintic 

(6) y 5 + a 4 ?/ 4 + a s y 3 + a 2 y 2 + a\y + x = 0. 

Computing the annihilating operator for the solutions of this equation has 
turned out to be a task of considerable computational complexity. The full 
output of the algorithm is a vector of five polynomials with 4306 monomials 
in total and is too large to display. The degrees of these polynomials 
with respect to the variables a\, . . . , 04, x are 15, 20, 21, 22, 23. The leading 
coefficient of the annihilating operator has degree 7 with respect to x and 
splits into the product of two factors. One of them is the discriminant 
of (ED while the other is a polynomial of total degree 15 with 264 terms. 
The largest of the numeric coefficients in the annihilating operator for the 
generic monic quintic equals 2739594525000. 

Example 11. A monic tetranomial with generic coefficients. By Theorem 
O the determination of a solution to the algebraic equation 

(7) y 6 + ay 2 + by + x = 

equals five. The roots of ([ID at a generic point x G C span the space of 
holomorphic solutions of the following fifth order linear differential operator 
with polynomial coefficients: 

(-255664128 a 10 + 395740000 a 5 6 4 + 1599609375 b 8 + 148780800 a 6 b 2 x + 2859609375 ab 6 x- 

499654656 a 7 x 2 - 1573425000 a 2 6 4 x 2 + 1051704000 a 3 b 2 x 3 + 16796160 a 4 x 4 ) 

(256 a 5 b 2 + 3125 b 6 - 1024 a 6 x - 22500 ab 4 x + 43200 a 2 b 2 x 2 - 13824 a 3 x 3 - 46656 x 5 ) £5+ 

(916300234752 a 16 + 18677130035200 a 11 6 4 - 38094525000000 a 6 b 8 - 134905517578125 ab 12 - 
77437887971328 a 12 b 2 x + 107910691200000 a 7 b e x + 332702753906250 a 2 b w x+ 
37877629059072 a 13 x 2 - 406052352000 a 8 b A x 2 + 552267618750000 a 3 b 8 x 2 - 
128020162314240 a 9 b 2 x 3 - 727448202000000 a 4 6 6 x 3 + 267464667561984 a 10 x 4 + 
43693344000000 a 5 b 4 x 4 - 1306049062500000 b 8 x 4 - 221398918963200 a 6 b 2 x 5 - 
2201395927500000 ab 6 x 5 + 359825022517248 a 7 x 6 + 1137850610400000 a 2 b 4 x e - 
711490376448000 a 3 b 2 x 7 - 10579162152960 a 4 x 8 )-£z+ 

30(-3264411795456 a 12 b 2 + 5653930240000 a 7 b % + 24016248046875 a 2 6 10 + 
3132022849536 a 13 x - 5896377011200 a 8 b 4 x - 145456875000 a 3 b 8 x- 
2678330105856 a 9 b 2 x 2 - 31473123300000 a 4 b 6 x 2 + 38750783864832 a 10 x 3 - 
41250841920000 a 5 6 4 x 3 - 221475515625000 6 8 x 3 - 21927996518400 a 6 6 2 x 4 - 
344319609375000 ab 6 x 4 + 52068310990848 a 7 x 5 + 163366013160000 a 2 6 4 x 5 - 
95796142732800 a 3 6 2 x 6 - 1311148560384a 4 x 7 )^ T 
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-120(-441539395584a 13 + 907546908800 a 8 b 4 + 3520234921875 a 3 6 8 - 
122881784832 a 9 b 2 x + 2644039462500 a 4 b 6 x - 14126136238080 a 10 x 2 + 
20651647680000 a 5 b 4 x 2 + 88754326171875 b 8 x 2 + 4397772787200 a 6 b 2 x 3 + 
122488079765625 ab 6 x 3 - 16700798121984 a 7 x 4 - 51772358925000 a 2 6 4 x 4 + 
28447624699200 a 3 b 2 x 5 + 352185242112 a 4 x 6 )^ 

-2520(4704614400 a 9 b 2 +4079375000 a 4 6 6 - 219268374528 a 10 x+ 
350501380000 a 5 b 4 x + 1464693750000 b 8 x + 5043513600 a 6 6 2 x 2 + 
1670568609375 ab 6 x 2 - 196541448192 a 7 x 3 - 597319515000 a 2 6 4 x 3 + 
304660958400 a 3 b 2 x 4 + 3302125056 a 4 x 5 )£+ 

5040(-2372960256 a 10 + 4421840000 a 5 b 4 + 18235546875 b s - 1113523200 a 6 b 2 x+ 
14120437500 ab 6 x - 1226244096 a 7 x 2 - 3581820000 a 2 b 4 x 2 + 
1587859200 a 3 6 2 x 3 + 13436928 a 4 x 4 ). 

Recall that the Newton polytope of a multivariate Laurent polynomial 
/Oi, . . . ,x n ) = 

a=(ai,...,a n )€AcZ n 



°a ^1 • • • 



supported in a finite set A is defined to be the convex hull of A. Intensive 
experiments suggest that there is an intrinsic relation between the two 
extreme coefficients in the optimal annihilating differential operator for 
an algebraic function. As we have seen in several examples before, the 
leading coefficient in the annihilating operator is typically given by the 
product of the discriminant of the defining algebraic equation and some 
other factor which has no apparent relation to the initial algebraic equation. 
The following conjecture suggests the structure of the Newton polytope of 
this polynomial factor. 



Conjecture 12. Let 



d d k 



dx k 

k=£ 

be the optimal annihilating operator for the algebraic function defined by 
the relation P(x, y) := y m + a\y mi + . . . + a n y mn + x = 0. Denote by 
) the discriminant of P(x, y) computed with respect to y. 
Then the polynomials 

Pd(ai, . . . ,a n ,x)/Q(ai, ...,a n ,x) and pe(ai, ...,a n ,x) 

consist of monomials with the same exponent vectors. In particular, they 
contain equally many monomials and have equal Newton polytopes. 



LINEAR DIFFERENTIAL OPERATORS FOR GENERIC ALGEBRAIC CURVES 13 

The following table summarizes the results of our computer experiments 
and illustrates Conjecture d21 It gives the linear ordinary differential op- 
erator whose solution space is spanned by the branches of an implicitly 
defined algebraic function y = y(x), the order of this operator and the 
multidegree of its leading coefficient with respect to x and the parameters 
of equation listed in lexicographic order. 
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Table 1: Computation times for annihilating operators and their properties 



Algebraic curve 


The annihilating operator for y = y(x) 


Order 


Leading 
coeff. 


(Jomput. 
time (sec.) 


y 4 + ay 3 + x = 


(27a 4 x 2 - 256x 3 )^ + 4x(27a 4 - 416x)^ 
+60(a 4 -36x)^- 360^| 


4 


(2,4) 


0.374 


y 4 + ay 3 + by 2 + x = 


x(14& a - Aa z b z + 8bx - 3a^x)(16& 4 - 4a^o 3 - 

1286 2 x + 144a 2 6x - 27a 4 x + 256x 2 )^ 
+ . . . + 120(746 3 - 21a 2 o 2 + 246x - 9a 2 x)^ 


4 


(4,6,7) 


1.03 


y 4 + ay 3 + by 2 + 
cy + x = 


(45c z + Ub 3 - 47a6c - 4a z b 2 + 12a 3 c + 8bx 
-3a 2 x)(27c 4 + . . . - 256c 3 ) 
. . . + 120(243c 2 + 746 3 - 249a6c - 21a 2 6 2 
+63a 3 c + 246x- 9a 2 x)^ 


4 


(4,7,7,6) 


3.011 


y 5 + ay + x = 


(256a 5 + 3125x 4 )^ + 31250x 3 ^+ 
73125x 2 ^ + 31875x^ - 1155 


4 


(4,5) 


2.012 


y 5 + ay 2 + oy + x = 


(512006° - 15930a 4 6 3 - 2187a 8 + 68000ao 4 x 
-1350a 5 6x - 26250a 2 o 2 x 2 )(256o 5 - 
27a 4 6 2 - 1600a6 3 x + 108a 5 x + 2250a 2 6x 2 

+3125x 4 )^ + . . . 
+120(4928006 6 - 139995a 4 6 3 - 16038a 8 + 
222000ao 4 x + 8100a 5 6x - 39375a 2 6 2 x 2 ) 


4 


(6,13,11) 


5.008 


y 5 + ay 3 + 6y 2 + cy 
+x = 


(102400c b + . . . - 2500a 3 6x 3 )(256c a - 276V 
+144a6 2 c 3 + . . . + 3125x 4 )-f^ 
+ . . . - 120(985600c 6 + . . . - 625a 3 6x 3 ) 


4 


(7,15, 
13,11) 


31.668 


y 5 + ay 4 + x = 


(265a 5 x 3 + 3125x 4 )£^ + 5x 2 
(384a 5 + 6875x)^r + 780x(4a 5 + 125x) 
+60(14a 5 + 1175x)^ + 6720^ 


5 


(4,5) 


10.436 


y 5 + ay 4 + fry 3 + x = 


(1680a6 y + . . . + 56a y x z )(108o i V- 
27a 2 6 4 x 2 + 22506 2 x 3 - 1600a 3 6x 3 + 

256a 5 x 3 + 3125x 4 )^ + ...+ 
1680(20160ao 9 + . . . + 22Aa 9 x 2 )± 


5 


(6,15,14) 


38.563 


y 5 + ay 4 + by 3 + 
cy 2 + x = 


(160380c 8 + . . . - 56a y x 3 )(108c 5 x+ 
166 3 c 3 x + ... + 3125x 4 )^ 
+ . . . + 1680(2779920c 8 - 1242720o 3 c 6 + . . . 
-224a 9 x 3 )^ 


5 


(7,16, 
14,13) 


279.242 


y 5 + ay 4 + 6y 3 + 
cy 2 + dy + x = 


(4928000d b + . . . + 56a y x 3 j(256d 5 - 27c 4 d z 

+1446c 2 d 3 + . . . + 3125x 4 )^§r 
+ . . . + 1680(901 1200d 6 + . . . + 224a 9 x 3 )^ 


5 


(7,17,15, 
13,11) 


3038.47 


y 6 + ay 3 + 6y 2 + 
cy + x = 


(6398437500c iU + . . . + 67184640a z 6c 2 x 5 j 
(3125c 6 + 25 66 5 c 2 + . . . 
-46656x 5 )J^ + ...+5040 
(72942187500c 10 + . . . + 13436928a 2 oc 2 x 5 ) 


5 


(10,18, 
17,16) 


799.427 
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